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ABSTRACT 

We analyse the magnetic activity characteristics of the planet hosting Sun-like star, HD 1237, using HARPS spectro-polarimetric 
time-series data. We find evidence of rotational modulation of the magnetic longitudinal field measurements consistent with our ZDI 
analysis, with a period of 7 days. We investigate the effect of customising the LSD mask to the line depths of the observed spectrum 
and find that it has a minimal effect on shape of the extracted Stokes V profile but does result in a small increase in the S/N (~ 

7%). We find that using a Milne-Eddington solution to describe the local line profile provides a better fit to the LSD profiles in this 
slowly rotating star, which also impacts the recovered ZDI field distribution. We also introduce a fit-stopping criterion based on the 
information content (entropy) of the ZDI maps solution set. The recovered magnetic field maps show a strong (4-90 G) ring-like 
azimuthal field distribution and a complex radial field dominating at mid latitudes (~45 degrees). Similar magnetic field maps are 
recovered from data acquired five months apart. Future work will investigate how this surface magnetic field distribution impacts the 
coronal magnetic field and extended environment around this planet-hosting star. 
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1. Introduction 

Observational studies of magnetism in late-type stars have 
evolved dramatically during the last two decades; from the clas¬ 
sical chromospheric activity diagnostics (e.g. Mount Wilson H- 
K project, Baliunas et al. 1995) to spectro-polarimetric snapshot 
surveys (the BCool project, Marsden et al. 2014) and detailed 
long-term magnetic monitoring (e.g. Morgenthaler et al. 2012). 
This has been enabled by the advent of improved instrumentation 
(e.g. ESPaDOnS@CFHT, Donati 2003) together with advanced 
data analysis techniques for detection (e.g. Least Squares Decon¬ 
volution, Donati et al. 1997) and mapping (e.g. Zeeman Doppler 
Imaging, Donati & Brown 1997) of magnetic fields. These recent 
studies have opened new possibilities for different areas of as- 
trophysical research, in particular, on dynamo processes and the 
origin of stellar magnetic fields across the HR diagram (see Do¬ 
nati & Landstreet 2009). In the case of main sequence solar-type 
stars, complex and relatively weak (« 1 kG) large-scale surface 
field topologies have been reported. The appearance of ring-like 
structures of significant (and even dominant) toroidal fields in 
these stars seems to be connected with the rotation period, and 
therefore with the dynamo mechanism behind their generation 
(Petit et al. 2008). 

Prom the perspective of exoplanet studies, activity-related 
signatures in late-type stars are known to affect the detection 
techniques in the form of radial-velocity jitter and photometric 
flicker (Bastien et al. 2014). By simulating the effects induced 
by active regions and spots in Sun-like stars, recent tools have 
been developed to estimate and remove their contribution from 
the observations (e.g. Dumusque et al. 2014). Other studies have 
considered the detectability of planets around active cool stars, 
by modelling the stellar activity from recovered starspot and 
magnetic field maps (Jeffers et al. 2014a; Donati et al. 2014). 


It is clear that knowing the characteristics of the stellar magnetic 
field is crucial for addressing the presence of exo-planets in a 
given system (e.g. the case of e Eridani, Anglada-Escude & But¬ 
ler 2012; Jeffers et al. 2014b). 

In addition, the stellar magnetic field dominates the environ¬ 
ment around late-type stars. This includes transient events such 
as flares and coronal mass ejections (Shibata & Magara 2011), 
and the development of persistent solar-like winds and astro- 
spheres (Wood 2004). These phenomena are known to have a 
profound impact on the structure of exoplanet atmospheres, a 
critical factor in the habitability of exo-planetary systems (see 
Vidotto et al. 2013; Cohen et al. 2011b, 2014). They can erode 
atmospheres through thermal evaporation, or non-thermal pro¬ 
cesses, such as sputtering and ion pick-up. Significant mass loss 
has been detected for exoplanets that is driven by the stellar 
wind (Vidal-Madjar et al. 2003; Linsky et al. 2010; Jensen et al. 
2012). Models of this interaction require accurate knowledge of 
the wind properties (Ekenback et al. 2010; Sanz-Eorcada et al. 
2011), and therefore, of the host-star surface magnetic field. 
However, given the observational limitations, robust surface field 
distributions are known for a very limited number of Sun-like 
planet-hosting stars (e.g. Eares et al. 2013; Eares 2014). In this 
context, detailed studies of these systems are very valuable re¬ 
sources not only as direct stellar counterparts of our solar system, 
but also for the growing interest in finding suitable Earth-like life 
supporting places in the Universe. 

In this article we present the detailed study of one such 
planet-hosting Sun-like star (HD 1237), in which we investi¬ 
gate the large scale magnetic field and chromospheric activity. 
This is the first step in characterising the impact the stellar mag¬ 
netic field on the circumstellar environment around this system. 
In particular, the conditions and possible interactions via the 
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magnetically-driven stellar wind, with the Jupiter-size exoplanet 
which comes as close as 0.25 AU in its orbit (Naef et al. 2001). 
In section 2, we summarise the main properties of the star. De¬ 
tails of the observations and calibration procedures are given in 
section 3. We present the activity diagnostics and variability in 
section 4. Section 5 contains a description of the implemented 
technique for extracting the magnetic field signatures from the 
spectro-polarimetric data. The required steps for the imaging 
procedure and the resulting surface field maps are presented in 
section 6. In section 7, we discuss our findings in context of pre¬ 
vious and ongoing studies of solar-type stars. Our main conclu¬ 
sions are summarised in section 8. 


2. HD 1237 stellar properties 

HD 1237 (GJ 3021) is abright (y„ag = 6.58), Sun-like star (G8V) 
located about 17.5 pc from the Sun, in the southern constella¬ 
tion of Hydrus (Koen et al. 2010). This object is a relatively 
young (~0.88 Gyr), chromospherically active, and confirmed 
exo-planet host star. Naef et al. (2001) used the enhanced chro¬ 
mospheric activity to explain the large residuals arising from 
the best Keplerian orbital solution of the planet (Mp sin(/) = 
3.37 + 0.14 Mop, Porb = 133.7 + 0.2 days, e = 0.51 + 0.02, 
a = 0.49 AU). 

Table 1 contains the basic stellar properties of HD 1237 taken 
from Ghezzi et al. (2010), Torres et al. (2006) and Saffe et al. 
(2005). Rotation period (Prot) estimates are sparse, ranging from 
~4.0 to 12.6 days, as summarized by (Watson et al. 2010). As 
presented in section 6.1, we obtain P^t = 7.0 + 0.7 from our 
observations (Sect. 3). In addition, we estimated the radial ve¬ 
locity, vr and the rotational velocity, v sin i, using an automatic 
spectral classification tool (MagIcS, Donati et al. 2012) and the 
fundamental properties of the star (Table 1). Using several of our 
observed spectra, we found on average a v sin i of 5.3 +1.0 km 
s“' and a vr of -5.2 + 0.2 km s"'. Literature values for vsin/ 
range between ~ 4.5 to 5.5 km s“^ (Naef et al. 2001; Torres 
et al. 2006; Schroder et al. 2009). However, as indicated by Naef 
et al. (2001 ), the v sin / value can be over-estimated for metal-rich 
stars such as HD 1237. For the subsequent analysis we adopted 
a V sin i of 5.3 km s“^ which is consistent, within the errors, with 
the value reported by Naef et al. (2001) and Torres et al. (2006). 


Table 1. HD 1237 basic properties. 


Parameter 

Value 

Reference 

S. Type 

G8V 

Torres et al. (2006) 

Teff [K] 

5572 + 40 

Ghezzi et al. (2010) 

log(g) 

4.58 + 0.2 

Ghezzi et al. (2010) 

R* [Rq] 

0.86 + 0.07 

Ghezzi et al. (2010) 

M. [Mq] 

1.0+ 0.1 

Ghezzi et al. (2010) 

vsini [km s“*]^ 

5.3 + 1.0 

This work 

Vr [km s“*] 

-5.2+ 0.2 

This work 

Prot [days] 

7.0 + 0.7 

This work 

log(Lx) 

29.02 + 0.06 

Kashyap et al. (2008) 

Age [Gyr]* 

~ 0.88 

Saffe et al. (2005) 


t: Other reports include 4.5 km s * (Schroder et al. 2009), 5.1 ± 
1.2 km s“‘ (Torres et al. 2006) and 5.5 ± 1.0 km s“’ (Naef et al. 
2001 ). 

L Age estimates range from 0.15 to 0.88 Gyr using various meth¬ 
ods (see Naef et al. 2001 and Saffe et al. 2005). 0.88 Gyr corre¬ 
sponds to the age determined using isochrones. 


3. Observational data 

We obtained observations using the polarimetric mode 
(Piskunov et al. 2011) of the HARPS echelle spectrograph 
(Mayor et al. 2003) at the ESO 3.6 m telescope at La Silla Obser¬ 
vatory. The wavelength coverage of the observations range from 
378 nm to 691 nm, with a 8 nm gap starting at 526 nm. 

Data were reduced using the REDUCE package (Piskunov 
& Valenti 2002; Makaganiuk et al. 2011), which was modi¬ 
fied for the HARPS instrument configuration. This package pro¬ 
duces an optimal extraction of the bias-subtracted spectra after 
flat-fielding corrections and cosmic ray removal have been car¬ 
ried out. The continuum level is determined by masking out the 
strongest, broadest features (e.g. the Balmer lines) and then fit¬ 
ting a smooth slowly varying function to the envelope of the en¬ 
tire spectrum. Spectra are obtained with resolutions varying from 
95 000 to 113 000, depending on the wavelength, with a median 
value of 106 000. Uncertainties are derived for each pixel assum¬ 
ing photon statistics. The star was observed at two epochs sepa¬ 
rated by 5 months (July and December) in 2012. A summary of 
the observations is presented in Table 2. 


Table 2. Journal of observations. The columns contain the date, the cor¬ 
responding Heliocentric Julian Date (HJD), the start time of the obser¬ 
vations in UT, the exposure times, and the Stokes I peak Signal-to-Noise 
ratio (S/N). The rotational phase (O) listed in the last column is calcu¬ 
lated using the rotation period derived in this work (Py^t = 7.0 d). 


Date 

(2012) 

HJD 

(2400000+) 

UT 

^exp 

[s] 

Stokes I 
Peak S/N 

Phase 

(®) 

First epoch 

Jul 15 56123.359 

08:04:36 

3600.0 

955 

0.000 

Jul 16 

56124.442 

10:03:43 

3600.0 

1214 

0.155 

Jul 17 

56125.399 

09:01:10 

3600.0 

841 

0.291 

Jul 18 

56126.361 

08:06:35 

3600.0 

877 

0.429 

Jul 19 

56127.440 

10:01:33 

3600.0 

753 

0.583 

Jul 20 

56128.356 

08:00:16 

3600.0 

1098 

0.714 

Jul 21 

56129.374 

08:25:27 

3600.0 

911 

0.859 

Jul 22 

56130.438 

09:47:58 

4800.0 

1092 

1.011 

Jul 23" 

56131.397 

08:32:33 

6680.0 

660 

1.148 

Jul 31 

56139.314 

06:49:52 

4800.0 

926 

2.280 

Aug 02 

56141.303 

06:32:21 

5000.0 

1040 

2.564 

Second epoch [20.29+ 

rotation cycles since Jul 15 2012] 


Dec 04 

56265.045 

00:35:27 

2800.0 

1205 

0.000 

Dec 05 

56266.045 

00:35:21 

2800.0 

964 

0.142 

Dec 06 

56267.044 

00:34:40 

2800.0 

1018 

0.285 

Dec 07 

56268.044 

00:33:59 

2800.0 

722 

0.428 


t: The listed values correspond to two spectro-polarimetric exposures 
merged in a single observation due to bad weather conditions. 


The exposure times listed correspond to one circularly polarised 
spectram (Stokes V) which results from combining four individ¬ 
ual sub-exposures using the ratio method. As is explained in Do¬ 
nati et al. (1997), the polarization signal is obtained by dividing 
spectra with perpendicular (orthogonal) polarization states (for 
HARPSpol and Stokes V; 45°, 135°, 225° and 315°, using the 
quarter waveplate). Additionally, a null-polarisation spectrum is 
constructed in order to check for possible spurious polarisation 
contributions in the observations. More details can be found in 
Bagnulo et al. (2009). Due to bad weather conditions, two con¬ 
secutive Stokes V spectra were added together for the night of 
2012 Jul23. 
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4. Magnetic activity and variabiiity 

In order to characterise the chromospheric activity level of the 
star during the observed epochs, we use the Call H (396.8492 
nm) & K (393.3682 nm) lines and the classic Mount Wilson S- 
index, S mw, defined as 


H and K represent the fluxes measured in each of the Ca line 
cores using 0.105 nm wide spectral windows. R and V are the 
fluxes measured in the continuum over 2 nm windows centred at 
390.1 nm and 400.1 nm respectively, on both sides of the Call 
region. 

4.1. Index calibration 

To compare the activity level of HD 1237 with other stars, we 
need to convert the measured HARPS S-index, S h, to the Mount 
Wilson scale. For this, we require a calibration factor, a, which is 
an instrument-dependent quantity that linearly relates the values 
for the classic S mw and the HARPS fluxes H, K, R, V: 



Sh 


We estimated a by including a set of stars with previous mea¬ 
surements of chromospheric activity via Smw (Santos et al. 
2000) and within the HARPS observations database. The spec¬ 
tral type of the reference stars and their reported S mw values are 
listed in the Table 3. The linear relation between Smw and the 
HARPS fluxes is plotted in the Figure 1 . The derived calibration 
factor, within the 1-cr uncertainty, is 

a = 15.39 + 0.65. (3) 

All spectra were co-aligned using a high S/N HARPS solar spec¬ 
trum as reference. We estimate a 5% typical error size in our 
HARPS flux measurements based on the possible differences in 
the continuum normalisation, which was performed in the same 
way for all the stars in the calibration. Therefore, the errors in 
S H are dominated by the conversion procedure. 

Table 3. Stars included in the a calibration. 


Name 

S. Type 

S MW 

C+MW 

HD 1835 

G3V 

0.364 

0.024 

HD 10700 

G8.5V 

0.173 

0.004 

HD22049 

K2Vk 

0.515 

0.026 

HD23249 

Kim-IV 

0.150 

0.012 

HD26965 

G9III-IV 

0.208 

0.018 

HD30495 

G1.5V 

0.292 

0.016 

HD61421 

E5IV-V 

0.187 

0.010 

HD76151 

G3V 

0.262 

0.018 

HD115617 

G7V 

0.161 

0.003 

HD 149661 

K2V 

0.356 

0.042 

HD152391 

G8.5Vk 

0.392 

0.030 

HD155885 

KIV 

0.400 

0.020 


We proceed with the estimation of the activity index S h, with 
the corresponding indicators Rhk (Middelkoop 1982) and 
(Noyes et al. 1984), which account for colour and photospheric 
correction respectively. 



0.005 0.010 0.015 0.020 0.025 0.030 

(H + K/R+V)h 


Fig. 1. Linear fit between Santos et al. (2000) Smw values and the 
HARPS fluxes A 5% error is estimated in our HARPS flux 

measurements. 

4.2. Activity indicators 

Figure 2 shows the cores of the Call H and K lines in the HARPS 
normalized spectra of HD 1237 for different observations, which 
are compared with the mean profile (red) derived from the entire 
data set. For the averaging procedure we take into account the 
slight differences in the wavelength range from each observation 
by an interpolation procedure to match the largest wavelength 
data points in the observed spectra. A high S/N HARPS Solar 
spectrum’, is shown in purple as reference. 

The upper plot shows one observation at a later epoch (2012 
Dec 07), where it is possible to observe a variation in the line 
profile. This is interpreted as a slight change in the chromo- 
spheric/photospheric activity of the star in comparison with the 
mean behaviour of the red line (that is dominated by profiles 
from the first epoch), especially in the H line (Fig. 2, right panel). 
The K-line region of the spectrum contained more noise. No 
change in the activity level of the star is visible in this particular 
line. 

Figure 3 shows the measured HARPS fluxes, (H+K)/(R + V), 
for each observation starting from 2012 Jul 15 (vertical blue line, 
HJD = 2456123.5). The x-axis units correspond to days after this 
initial date. The activity of the star showed a marginal variation, 
in a similar way as in the first epoch. This may be due to the 
rotation of active regions over the stellar surface. The similarities 
between the activity levels between both epochs, could be an 
indication of a stable large-scale magnetic field configuration. 

The chromospheric activity level of the star can be quanti¬ 
fied by using Eq. (2), and the derived value for a in Eq. (3). We 
estimate an average value of 5 h - 0.46 + 0.02 for the available 
observations of HD 1237. Similar activity levels have been re¬ 
ported for the Sun-like star ^ boo A (Tefr = 5600 K, Prot = 6.4 
days. Age: ~ 0.2 Gyr, Mamajek & Hillenbrand 2008; Morgen- 
thaler et al. 2012). As a reference value, the Solar S-index is 


' S/N: 347 @550 nm. Date: 2007 Apr 12 - Low activity period. 
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X [nm] 
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Fig. 2. Core regions of the Call K (left) and H (right) lines of HD 1237. Four spectra from our sample are plotted, with the corresponding dates on 
the left y-axis and vertically shifted (0.25 units) for visualisation purposes. The black spectrum at the top was taken at a much later epoch (2012 
Dec 07). The red line shows the mean profile for the entire dataset, while the purple line is a HARPS solar spectrum used as reference. 


Sq ^ 0.1783, with a variation of ~ 0.02 from solar maximum to 
minimum (Lockwood et al. 2007). 

Using the mean derived values of the S - index we can now 
apply a transformation to obtain the parameter /?hk. which takes 
into account the colour of the star in the activity estimation (Mid- 
delkoop 1982). /?hk is defined as 

^hk = (Ccf)(5h)(1.34x10-4), (4) 

where Ccf is a colour-dependent function. For main sequence 
stars with 0.3 < (B -V) < 1.6, Cqf is given by 

log(CcF) = 0.25(B-y)^-1.33(B-y)2-i-0.43(B-y)-i-0.24 . (5) 

Noyes et al. (1984) derived an expression in order to include 
photospheric corrections to the Call core fluxes, which is 
written as 

^FlK ~ ^FIK ~ ^phot . (6) 

with f^phot expressed also as a function of (B - y): 

log(/?phot) = -4.898 -H 1917(B - vf - 2.893(B - y)^, (7) 

and valid in the range of 0.44 < (B - y) < 1.0. Table 4 sum¬ 
marises the activity indicators for the available observations of 
HD 1237. The errors quoted are the mean measurement errors 
for these quantities. B and V magnitudes were taken from the 
Hipparcos catalog (B = 7.335 mag, V = 6.578 mag, Koen et al. 
2010). High activity levels are related with magnetic fields on the 
stellar surface. These magnetic signatures are encoded in the po¬ 
larised spectra of the star and will be discussed in the following 
section. 


Table 4. Average activity indicators for HD 1237. 


<5h) 

log(RHK) 

log(^HK) 

0.462 + 0.019 

-4.29 + 0.04 

-4.38 + 0.05 


5. Magnetic field signatures 

The signal-to-noise ratio (S/N) in the observations is not suffi¬ 
ciently high to detect magnetically-induced spectro-polarimetric 
signatures in single lines. However, by applying a multi-line 
technique, e.g. Least Squares Deconvolution (LSD, Donati et al. 
1997), it is possible to increase the S/N by a factor of ~50 - 100, 
adding-up the signal from thousands of spectral lines over the 
entire spectral range (for HARPS: 378-691 nm), enhancing 
our sensitivity to magnetic signatures in the observations (see 
Kochukhov et al. (2010) for a recent review of the LSD tech¬ 
nique). This procedure requires a photospheric model (line list) 
matching the spectral type of our target star. This is done using 
an atomic line list database ’ (Kupka et al. 2000) and the stellar 
fundamental properties listed in table 1. We assumed a micro¬ 
turbulence parameter of 1.3 km s“' (Ghezzi et al. 2010) and solar 
abundance for the photospheric line list which included ~ 15000 
lines within the HARPS spectral range. 

From this initial photospheric line list we generated two dif¬ 
ferent masks used with the LSD calculation. In the first mask the 
strong lines (and lines blended with these lines) that form in the 
chromosphere or that break the basic assumptions of LSD, (e.g., 
Ca II H&K, Ha) are removed (cleaned mask). This reduced the 
number of lines included to ~ 11000. After the mask cleaning, 
a numerical routine based on the Levenberg-Marquardt, non¬ 
linear least-squares algorithm from the mpfit library (More 1978; 
Markwardt 2009) is applied to fit the line mask to the observed 

"http://vald.astro.uu.se/ - Vienna Atomic Line Database 
(VALD3) 


Article number, page 4 of 13 


















































Alvarado-Gomez et al.: Activity and Magnetic Field Structure of HD 1237 


0.034 


0.032 


> 

+ 

cr 

> 0.030 

+ 

I 


0.028 


0.026 


Fig. 3. HARPS fluxes, (H + K)/(R + V), for the available observations of HD 1237. The blue and green vertical lines denote the beginning of each 
observed epoch. The x-axis contains the number of days since the Heliocentric Julian Date (HID = 2456123.5) of the first observation (2012 Jul 
15). The last two data points of the first epoch are unevenly distributed (after the vertical black line). We estimate a 5% error in our measurements. 
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Fig. 4. Calculated Bi for the available observations of HD 1237. The blue and green vertical lines denote the beginning of each observed epoch. 
The x-axis contains the number of days since the Heliocentric Julian Date (HJD = 2456123.5) of the first observation (2012 Jul 15). The last two 
data points of the first epoch are unevenly distributed (after the black line). Colours indicate the line mask used for the LSD procedure. 


Stokes I spectrum, adjusting the individual depths of the spectral 
lines (cleaned-tweaked mask). This is performed through the en¬ 
tire HARPS wavelength coverage. This step is more commonly 
carried out when applying LSD to hot (OB-type) stars which 
have fewer lines (e.g, Neiner et al. 2012). Finally, LSD was ap¬ 
plied to the spectro-polarimetric data using the final masks, gen¬ 
erating in this way a single, averaged line profile per observa¬ 


tion (LSD Stokes I, V and diagnostic N profiles). A velocity step 
Av = 1.4 km s ' was used to construct the LSD profiles. This 
velocity spacing considers two pixels per spectral element of the 
instrument (in the case of HARPSpol, R - 2.5 km s ' and 3.4 
px per resolution element). 

Both procedures lead us to consistent results in the obtained 
LSD signatures of the star. While no clear change is observed in 
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the null polarisation check, subtle qualitative differences appear 
in both Stokes profiles, in the sense that the cleaned-tweaked 
mask seems to get a broader unpolarised profile with a slightly 
weaker signature in the circular polarised profile, in comparison 
with the clean line mask. These minimal shape differences in 
the LSD line profiles can have a much larger effect in hot stars 
where fewer lines are generally available. On the other hand, the 
resulting S/N of the LSD profiles for each individual observation 
was systematically higher using the cleaned-tweaked mask than 
in the cleaned case, despite the same number of spectral lines in 
their masks (11048). On average, a ~7% increase was obtained 
in the S/N of the LSD profiles with the cleaned-tweaked mask. 

5.1. Longitudinal magnetic field 

Using the derived Stokes 1 and V LSD profiles, it is possible to 
obtain information about the surface averaged longitudinal mag¬ 
netic field {B(). Following Donati & Landstreet (2009), an esti¬ 
mate of B( (in G) is given by 


B{ = -714 


f vN(v)dv 
Ag][\-l{v)]dv ’ 


( 8 ) 


where the radial velocity shift v (in km s ') is measured with 
respect to the average line derived from LSD, with central wave¬ 
length A (in pm) and mean Lande factor g. As this measurement 
is an integrated quantity over the visible surface, it cannot pro¬ 
vide complete information for stars that host complex large-scale 
magnetic fields. From multiple measurements of Be taken over 
a stellar rotation period it is possible to gain a first insight con¬ 
cerning inhomogeneities of the disk-integrated magnetic field. 
Hence, it is also possible to estimate the stellar rotation period, 
using the modulation in a time-series of Be measurements, pro¬ 
vided that they span more than one rotation period. 

Figure 4 shows the measurements of Be for the available ob¬ 
servations of HD 1237. The integration was centered in the radial 
velocity of the star (-5.2 km s“'), and covering the entire Stokes 
V signature (+12.5 km s“*). The uncertainties were computed 
from standard error propagation from the spectra. The colours 
denote the line mask used in the LSD procedure, showing that 
the slight differences in each of the measurements are consistent 
within the errors. This implies that either set of LSD profiles can 
be used to obtain robust longitudinal magnetic field measure¬ 
ments. As the two sets of profiles are so similar, there is only 
a difference in S/N without any noticeable impact on the maps 
structure. 

Be shows a clearer rotational variation in comparison with 
S H (Fig. 3), reflecting the additive nature of the chromospheric 
emission, in contrast with the mixed polarity effects of the mag¬ 
netic field. The behaviour of the longitudinal field is consis¬ 
tent through the observed epochs, with an varying amplitude of 
roughly ~10 G. This value is somewhat larger in comparison 
with the solar value (Be < 4 G, Daou et al. 2006; Kotov et al. 
1998) and the average value from snapshot observations of other 
Sun-like stars of the same and different spectral types (3.3, 3.2 
and 5.7 G for F, G and K-dwarfs respectively, Marsden et al. 
2014). However, caution is advised in these averaged compar¬ 
isons given the rotational variability of Be and the nature of a 
snapshot survey. Similar variations have been observed in the 
long-term monitoring of the active Sun-like star ^ Boo A (Mor- 
genthaler et al. 2012). 


6. Surface magnetic field mapping 

6.1. Optimal line profile and stellar parameters 

We reconstructed surface magnetic field maps by applying the 
tomographic inversion technique of Zeeman Doppler Imaging 
(ZDI) (Vogt et al. 1987; Semel 1989). ZDI, described by Hus¬ 
sain et al. (2000), has been used to recover magnetic field maps 
on the surfaces of stars, ranging from T Tauri stars to binary sys¬ 
tems (e.g. Barnes et al. 2004; Dunstone et al. 2008). The code 
recovers the magnetic flux distribution across the stellar disk, 
modulated by the stellar rotation, by using time series of pho- 
tospheric absorption line profiles (LSD Stokes I) and circularly 
polarised profiles (LSD Stokes V). 

In order to recover reliable magnetic field maps, it is neces¬ 
sary to properly model the LSD Stokes I and V profiles and their 
temporal variation. To do this, good constraints should be ob¬ 
tained on the local line profile description and the stellar param¬ 
eters. Two different synthetic line shapes were tested; a Gaus¬ 
sian profile, which is commonly used in magnetic field studies 
in Sun-like stars (e.g. Boro Saikia et al. 2015), and a Milne- 
Eddington profile, fitted to a solar LSD profile derived from a 
high S/N HARPS spectrum. This last approach was previously 
considered in ZDI of accreting T Tauri stars (Donati et al. 2008a) 
and M-dwarfs (Morin et al. 2008). For both cases, we assumed 
a linear dependance of the continuum limb darkening with the 
cosine of the limb angle (slope u ^ 0.65, Sing 2010). 

In addition, we estimated the rotational period Prot, differen¬ 
tial rotation profile (e.g. Q(/) = Qeq “ t/Gsin^(/), see Petit et al. 
2002), and inclination angle i of the star. This is done by gen¬ 
erating a grid of ZDI models, covering a range of values for the 
involved quantities, and minimising the reduced from syn¬ 
thetic line profile fitting (see Collier Cameron 1995; Hussain 
et al. 2009). Figure 5 (left) shows the results of the minimisa¬ 
tion analysis over Pjof As was mentioned in Sect. 2, the rotation 
period of HD 1237 is not well known. Our analysis shows the 
harmonic behaviour of this parameter, with a fundamental value 
of ~ 7.0 days. We estimate a 10% error for this period determi¬ 
nation, given the width of the minima in the left panel of Fig. 4. 
This estimate is consistent with the high activity level measured 
in this star (Sect. 4). For completeness, we performed a Lomb- 
Scargle periodogram, using the first epoch dataset, on the tempo¬ 
ral variations of Be (Fig. 4). Figure 5 (right) shows the obtained 
power spectrum, having a best-fit period of Pjot = 6.8 ± 0.2 days 
with an associated P-Value statistic of 5.18 x 10“^ (Zechmeister 
& Kiirster 2009). 

No good constraints were obtained by x^ minimisations 
for the differential rotation parameters and the inclination an¬ 
gle of the star. Therefore, no differential rotation profile was 
included in the mapping procedure. For the inclination angle, 
we considered the expected value from rigid body rotation, i.e. 
sin/ = (Prot ■ vsin/)/(27rPt). Given the uncertainties of the in¬ 
volved quantities (see Table 1), the inclination angle should lie 
somewhere between ~ 40 - 60°. We performed reconstructions 
for 40, 50 and 60 ° and find no substantial differences between 
these magnetic field reconstructions. We present here the maps 
obtained assuming a 50 ° inclination angle. 

Finally, we compared the optimum stellar and line parame¬ 
ters associated with the LSD profiles produced using both masks 
(cleaned mask and clean-tweaked mask). No significant differ¬ 
ence is found in the optimal parameters. Hence in the subsequent 
analysis we only consider the LSD profiles produced using the 
cleaned-tweaked mask. 
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Fig. 5. Left: Minimisation results for HD 1237 showing the reduced as a function of Prot- The black segmented lines corresponds to a fourth- 
order polynomial fit. Right: Power spectrum obtained from the Lomb-Scargle periodogram analysis on the temporal variation of Bf (Fig. 4). 


6.2. Optimal Fit Quality: Entropy Content in ZDI Maps 

The last of part of the analysis corresponds to the selection of 
the optimal fit quality (reduced of the model with respect 
to the observations. In principle, the recovered field distribution 
and associated profiles should try to achieve the lowest possi¬ 
ble value of x^- Still, given the limitations of the observations 
and the ZDI technique (S/N, spatial resolution, phase-coverage, 
etc.), the goodness-of-fit level has to be determined carefully to 
avoid the appearance of numerical artifacts in the final maps. 
This is particularly important in the case of resolution-limited 
maps (i.e. slowly rotating stars). However, there are few pub¬ 
lished procedures to establish a robust “stopping criterion”, i.e. 
the point at which noise starts to affect the reconstructed image. 
Motivated by this, we propose a systematic method for estimat¬ 
ing the optimal fit quality for a given set of ZDI magnetic field 
maps, using HD 1237 as a test-case. It is important to note here 
that since this procedure is defined a posteriori over the resulting 
maps themselves (2D images), it does not modify the regulari- 
sation functions imposed to ZDI. In this sense, its application to 
other stellar systems should be straightforward. 

We begin with the definition of the entropy 5, as an estimate 
of the information content in an image*. Following Sonka et al. 
(2007), let P{k) be the probability that the difference between 
two adjacent pixels is equal to k. The image entropy can be esti¬ 
mated as 

5 =-J]p,log2(P,), (9) 

k 

where log 2 is the base 2 logarithm. This implies that a larger or 
smaller amount of entropy in the image will depend on the con¬ 
trast between adjacent pixels. An image that is perfectly constant 

^ A similar implementation of entropy is commonly used as a regu- 
larisation function in ZDI (see Piskunov & Kochukhov 2002). 


will have an entropy of zero. For the methodology described be¬ 
low, we are not interested in the absolute values of the entropy, 
rather its overall behaviour as a function of the reduced 

For a given converged ZDI solution (i.e. a particular value of 
reduced ;y^), we can calculate the total entropy content {S t) by 
applying the definition given by equation (9) to each one of the 
recovered maps: 

= + + ( 10 ) 

where Sr, Sm, and S a represent the entropy contained in the ra¬ 
dial, meridional, and azimuthal field components maps, respec¬ 
tively. Figure 6 shows in red the behaviour of 5 t as a function of 
the reduced for the ZDI solution set of HD 1237. Each point 
in Fig. 6 corresponds to a converged ZDI solution. 

As expected, by decreasing the reduced x^ the information 
content in the resulting ZDI solution increases (field strength 
and structure). The green and blue symbols indicate the entropy 
growth (first derivative) and its variation (second derivative), re¬ 
spectively. The overall behaviour of the total entropy content and 
information growth is consistent for both cases. In the case of the 
Gaussian profile (Fig. 6, left panel), the entropy growth remains 
fairly constant (close to zero) for large reduced values, reach¬ 
ing a maximum* around ^ 1.2. However, by that point the 
concavity of the curve has changed (negative second derivative), 
suggesting a different regime for the information growth in the 
ZDI solutions. This is interpreted as a noise signature, reflected 
as artifacts in the final maps leading to an additional increment 
in the information growth. For this reason we adopt as optimal 
fit level, the reduced x^ value for which the rate of change in the 
information growth (second derivative) in the ZDI solution set is 
maximised. In this particular case, this occurs around x^ - ^ -4, 

^The apparent negative sign in the first derivative is due to the re¬ 
versed direction of the x-axis (reduced ^^). 
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Fig. 6. Optimal selection criteria applied to the ZDI solution set of HD 1237. Each panel contains the results for the July dataset using the 
Gaussian (left) and Milne-Eddington (right) line profiles. The red symbols show the behaviour of the total entropy content 5 t as a function of the 
reduced Each point corresponds to a converged ZDI solution. Green and blue symbols represent the first and second derivatives as indicated. 
A fourth-order polynomial fit has been applied (segmented line) to find the value for which the rate of change in the information growth (second 
derivative) in the ZDI solution set is maximised. The optimal fit level is indicated in each case. 


as is indicated by the fourth-order polynomial fit (segmented 
line) in the left panel of Fig. 6. A similar criterion plot is con¬ 
structed for the Milne-Eddington line profile (Fig. 6, right panel), 
where a lower optimal reduced ^ 1.1 is achieved in this case*. 

6.3. ZDI Maps and Synthetic Stokes V Profiles 

We reconstruct the ZDI surface magnetic field maps and the syn¬ 
thetic circularly polarised profiles based on the time-series of 
LSD Stokes V spectra. For this we use the cleaned-tweaked line 
mask (Sect. 5), and the stellar and line parameters derived in 
Sect. 6.1 . The goodness-of-fit level in each case is selected un¬ 
der the criterion described in the last section. 

Figures 7 and 8 show the results of the ZDI procedure for 
the first epoch observations (2012 July) of HD 1237. Four verti¬ 
cal panels are presented, where the first three correspond to the 
Mercator-projected magnetic field maps in Gauss (G), decom¬ 
posed into the radial (Br - top), meridional (Bm - middle), and 
azimuthal (Ba - bottom) components. The phase coverage is in¬ 
dicated by the black tick marks in the upper y-axis. The spatial 
resolution of the maps is ~ 18° in longitude at the stellar equator 
(at the poles, the map has a considerably poorer resolution than 
the equator). However, the phase coverage also has a significant 
impact on the level of detail that can be recovered. In this case, 
the right hand side would have a slightly better resolution than 
the left hand side of the image. The last panel shows the fitted 
synthetic Stokes V profiles to the spectro-polarimetric observa¬ 
tions in each rotational phase (O). 

*The same criterion was applied to generate the magnetic field maps 
in the December dataset (Appendix A). An optimal reduced - 0.6 
was obtained in this case, as is expected for a dataset with fewer con¬ 
straints. 


The recovered maps show a relatively complex field distribu¬ 
tion across the surface. The field is dominated by the azimuthal 
component, displaying a strong (—h90 G) large-scale ring-like 
structure around 45 ° in latitude. In the radial component, two 
large and moderately strong (~ +50 G) regions of opposite polar¬ 
ities are also located at higher latitudes while weaker (~ +25 G) 
small-scale features of mixed polarities appear close to the equa¬ 
tor. These large magnetic features are somewhat preserved, with 
reversed polarities, in the meridional magnetic field maps. How¬ 
ever, some cross-talk from the radial component may be present 
in the meridional map (see Donati & Brown 1997). 

Although the overall large-scale structure and field strength 
are consistent between both line profiles, several differences are 
clearly visible in the characteristics of the field components. 
First, the maps recovered using the Milne-Eddington line shape 
include additional small-scale features leading to a more com¬ 
plex field distribution in the surface. These are more promi¬ 
nent in the maps for the radial and meridional components. 
Second, the recovered magnetic field distribution, in the Milne- 
Eddington case, seems to be slightly shifted to lower latitudes. 
These differences can be understood from the fact that the Milne- 
Eddington line shape is a better representation of the derived 
LSD Stokes I profile. This is also true for the the shape of the cir¬ 
cularly polarised profile, leading to a more detailed (additional 
small-scale structures) and a field distribution shifted to lower 
latitudes (as a consequence of the sensitivity in the core of the 
profile). Both elements are translated into the lower optimal re¬ 
duced value that can be reliably achieved in this case (see 
Sect. 6.2), and therefore an improved fit (lower panel in Fig. 8). 

The reconstructed map for the second epoch observations 
(2012 December) and the corresponding synthetic Stokes V pro¬ 
files are presented in the Appendix A. 
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Fig. 7. Results of the ZDI analysis for the first epoch observations of 
HD 1237 using the Gaussian line profile. The first three panels show 
the surface magnetic field components Br, Bm, and Ba, respectively. 
The colour scale indicates the polarity and the magnitude of the mag¬ 
netic field component in G, while the phase coverage is indicated by 
the black tick marks in the upper y-axis. The segmented horizontal line 
indicates the surface visibility limit, imposed by the adopted inclination 
angle of the star (i = 50 °). The last panel show the comparison between 
synthetic (red) and observed (black) Stokes V profiles obtained for this 
particular epoch, in each observational phase O, where the recovered 
maps fit the spectro-polarimetric data to an optimal reduced = 1.4. 


7. Summary and discussion 

In this paper, we presented a detailed study, covering two ob¬ 
servational epochs, of the activity and magnetic field structure 
of the Sun-like planet-hosting star HD 1237. The chromospheric 
activity level of HD 1237 estimated from the calibrated 5 h = 


Fig. 8. Results of the ZDI analysis for the first epoch observations of 
HD 1237 using the Milne-Eddington line profile. See caption of Figure 
7 for more details. In this case, the maps fit the spectro-polarimetric data 
to an optimal reduced = 1.1. 


0.46+0.02 and log(R'^^) — -4.38+0.05 values is similar to other 
active Sun-like stars (e.g. f boo A, Morgenthaler et al. 2012), 
and considerably higher than the Solar case (~ 3 times higher 
in terms of the average S-index, Lockwood et al. 2007). A much 
larger difference has been reported for the X-ray activity level, 
with a log(Lx) = 29.02+0.06 for HD 1237 (Kashyap et al. 2008), 
which is two orders of magnitude higher than the value estimated 
for the Sun during solar maximum (Peres et al. 2000). The chro¬ 
mospheric activity level remained fairly constant over the ~ 5 
month period between the two sets of observations. Our estimate 
falls between the previous measurements of log(/?)^^) = -4.27 
(Naef et al. 2001) and log(7;;^^) = -4.44 (Saffe et al. 2005). 
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Given the large uncertainties in the activity indicators, it is diffi¬ 
cult to address if these variations have some correspondence to 
a magnetic cycle in the star. 

In order to extract magnetic field signatures from spectro- 
polarimetric data, we applied the LSD multi-line technique to the 
observations. Two different line-masks for the LSD profiles were 
compared and tested. The standard procedure involves employ¬ 
ing a mask “cleaned” of chromospheric and strong NLTE line 
profiles (e.g. Marsden et al. 2014); we compared the results from 
this procedure with improving the line-list further by “tweaking” 
line-strengths, so that they are tailored to the line depths in the 
observed spectrum. Both approaches lead to similar results in 
the obtained LSD profiles, and therefore in the physical quanti¬ 
ties inferred (e.g. the longitudinal magnetic field, B(). However, 
for the same number of spectral lines, the average S/N of the 
LSD profiles recovered with the aid of the cleaned-tweaked line 
mask was ~ 7% higher than in the cleaned case. 

The longitudinal magnetic field {B(), estimated from the 
Stokes I and V LSD profiles, showed a clear rotational modula¬ 
tion with an amplitude of ~10 G. This behaviour was preserved 
in both observed epochs. Placing these measurements in con¬ 
text of other stars, similar variations have been observed in the 
long-term monitoring of the active Sun-like star f Boo A (~ 4-9 
G, Age; ~ 0.2 Gyr, Mamajek & Hillenbrand 2008; Morgenthaler 
et al. 2012), the K-dwarf exoplanet host e Eri (~ 10-12 G, Age: 
~0.2 - 0.8 Gyr, Janson et al. 2008; Jeffers et al. 2014b), and 
more recently for the young solar analogue HN Peg (~ 14 G, 
Age: ~0.2 Gyr, Eisenbeiss et al. 2013; Boro Saikia et al. 2015). 
The chromospheric and X-ray activity levels of these stars are 
also very similar in comparison to HD 1237. ^ Boo A has an av¬ 
erage S-index of 5 hk - 0.45 and log(Lx) - 28.91 (Gray et al. 
1996; Wood & Linsky 2010). e Eri has strong magnetic activ¬ 
ity with a mean 5 hk - 0.50 and log(Lx) - 28.22 (Jeffers et al. 
2014b; Poppenhaeger et al. 2011). Similarly, previous reports for 
HN Peg show a 5 hk - 0.35 and log(Lx) - 29.19 (Boro Saikia 
et al. 2015; Schmitt & Liefke 2004). Despite the various simi¬ 
larities among these systems, the complete relation between the 
magnetic field and its influence over different layers of the stellar 
atmosphere (activity) is not yet fully understood. Similar to these 
other systems, the chromospheric activity of HD 1237 does not 
show a clear correlation with Bf. This is interpreted as the result 
of probing different spatial and temporal scales in each of these 
measured quantities. 

We recovered the optimal stellar parameters of HD 1237 
using a reduced minimisation scheme, based on the tomo¬ 
graphic inversion technique Zeeman Doppler Imaging (ZDl). 
The analysis yields a rotation period, Prot = 7.0 + 0.7 days. Liter¬ 
ature values of Prot are uncertain, ranging between ~ 4.1 - 12.6 
days (see Watson et al. 2010 and references therein). The 7.0 
days value found in this work is consistent with the ZDI anal¬ 
ysis and both, the chromospheric and coronal activity levels of 
the star. Similar procedures were applied to estimate the inclina¬ 
tion angle and the differential rotation of HD 1237. However, the 
available observations did not provide enough constraints for a 
robust determination of these parameters. Therefore, we consid¬ 
ered an inclination angle derived using the stellar properties of 
the star and solid body rotation (i.e. i ~ 50 °). No differential ro¬ 
tation profile was included in the reconstruction of the magnetic 
field maps. 

Eor the surface magnetic field mapping procedure, two dif¬ 
ferent synthetic line shapes (Gaussian/Milne-Eddington) were 
tested in order to investigate their impact on the ZDI maps 
for slowly rotating solar-type stars. We showed that, both pro¬ 
files recover robust magnetic field maps. However, the Milne- 


Eddington line profile yields a better spectro-polarimetric fit 
(lower optimal reduced x^) leading to a more detailed structure 
recovered in the ZDI maps compared to the Gaussian case. In 
connection to this, a fit-stopping criterion based on the informa¬ 
tion content (entropy) of the ZDI maps solution set was intro¬ 
duced. This allows the identification of the optimal reduced x^ 
value, avoiding to some extent possible the appearance of numer¬ 
ical artifacts in the ZDI maps. The optimal fit level is given by 
the reduced x^ value for which the rate of change of the entropy 
growth in the ZDI solution set is maximised. 

The large-scale magnetic field of HD 1237 showed a com¬ 
plex distribution at the stellar surface. The strongest magnetic 
field features appear at mid latitudes (~ 45 °) in the azimuthal and 
radial components. The field is dominated by the azimuthal com¬ 
ponent, displaying latitudinal belts or ring-like structures across 
the stellar surface. This has been observed in other studies of 
active Sun-like stars (e.g. Eolsom et al. 2013) and in numerical 
simulations of a rapidly rotating Sun (Brown et al. 2010). As has 
been suggested previously, the appearance of significant, or even 
dominant toroidal fields in the surface of these types of stars is 
connected with the rotation period, with ~ 12 days as the rota¬ 
tion threshold (Petit et al. 2008). In comparison with other cool 
stars, HD 1237 {P^ot = 7.0 days, this work, M, = 1.0 + 0.1 Mq, 
Ghezzi et al. 2010) would be located near ^ Boo A (G8V, Age: 
~ 0.2 Gyr) and GJ 182 (MOV, Age: 0.50+^ ° Gyr, Liu et al. 2004) 
in the mass-period, large-scale magnetic field diagram (Eig. 3 in 
Donati & Landstreet 2009), with a mostly toroidal field topology 
and a considerable contribution from the radial component. With 
the same spectral type as HD 1237, ^ Boo A displays a slightly 
weaker field in its surface (~ +60 G), with an alternating dom¬ 
inance between the radial and the azimuthal components in the 
observed long-trend evolution (~ 4 years, Morgenthaler et al. 
2012). During the observed azimuthal-dominated epochs of this 
star, the field distribution is highly similar to the one derived for 
HD 1237 in this work. A strong uni-directional azimuthal field 
appears at low latitudes, with large mixed polarity regions in the 
radial field and a minor contribution from the meridional compo¬ 
nent. However, the magnetic regions of ^ Boo A are much larger 
and less concentrated as the ones found for HD 1237. This could 
be related with the small difference in their rotation periods (6.43 
days for ^ Boo A and 7.0 days for HD 1237) and/or with the dif¬ 
ferential rotation that may be occurring in the surface. A similar 
situation appears in the case of GJ 182 in terms of the recov¬ 
ered magnetic field topology. However, the field strength is much 
larger in this case (up to 400 G, Donati et al. 2008b), which is 
mostly connected with the rapid rotation, large latitudinal shear, 
and the nearly fully convective nature of this star. 

Regardless of the low vsin/ of HD 1237, it is clear that a 
more complex field distribution is required to fit the observed 
spectra. Some weaker small-scale regions are recovered closer 
to the stellar equator. Despite their relative low strength com¬ 
pared to the total surface field, the contribution from these small- 
scale features to magnetically-related phenomena may be signifi¬ 
cantly higher. These features, in combination with strong mixed- 
polarity regions missed by ZDI, can influence the quiescent coro¬ 
nal emission (Johnstone et al. 2010), the X-ray modulation (Ar¬ 
zoumanian et al. 2011), and the wind structure around, very ac¬ 
tive stars (Lang et al. 2014), and the Sun (Garraffo et al. 2013). 

The large-scale held distribution appears similar at both of 
the observed epochs, confirming the observed chromospheric ac¬ 
tivity and longitudinal magnetic held behaviour. Using the final 
ZDI maps, and the derived synthetic Stokes V profiles, we were 
able to reproduce the variations of Be consistently (for both ob¬ 
served epochs). As an example, the maximum value of Be in Eig. 
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4, coincides when the large positive polarity region of the radial 
field is located close to the limb (day 3). Similarly, the minimum 
value obtained for B[ in Fig. 4, results from the negative polar¬ 
ity region in the radial field map, this time located much closer 
to the disk centre. This corroborates the robustness of our ZDI 
magnetic maps for this system. 

As was mentioned in Sect. 2, HD 1237 has a planetary com¬ 
panion with a mean separation of 0.49 AU, and a projected mass 
of Mp sin(0 = 3.37 + 0.14 Mtj,. This exoplanet is relatively far 
out in comparison with other similar systems where ZDI maps 
of the host star are available (see Fares et al. 2013 and references 
therein). However, the enhanced activity levels of HD 1237 and 
the relatively strong and complex surface magnetic field (in con¬ 
nection with its stellar wind), could impact significantly the con¬ 
ditions experienced by the exoplanet through its orbit. Previ¬ 
ous parametric studies of this system have predicted a relatively 
large mass loss rate (~ 85 Mq), which could even lead to mag- 
netospheric radio emission from the exoplanet (Stevens 2005). 
This will be considered in a future study, incorporating the re¬ 
covered ZDI magnetic field maps into a detailed 3D magneto¬ 
hydrodynamics (MHD) code (BATS-R-US, Powell et al. 1999; 
Toth et al. 2012), originally developed and validated for the solar 
wind and corona (e.g. Sokolov et al. 2013; van der Holst et al. 
2014), and recently applied in the stellar context (e.g. Cohen & 
Drake 2014; Cohen et al. 2014). 

8. Conclusions 

We have reconstructed magnetic field maps of the young planet¬ 
hosting G-type star, HD 1237 using the technique of Zeeman 
Doppler Imaging. As part of this detailed spectro-polarimetric 
study we have found the following; 

o We find that assuming a Milne-Eddington approximation for 
the local line profile produces a better fit to the shape of the 
observed LSD profiles. This influences the amount magnetic 
structures that can be recovered in the magnetic field maps. 

o We propose a robust method to define the “stopping crite¬ 
rion” in Zeeman Doppler Imaging techniques. This allows 
one to choose the optimum degree of fit, beyond which the 
model ceases to provide a good fit to the observed dataset 
due to the introduction of artifacts into the resulting image. 
We have successfully applied this to two different datasets, 
with vastly different degrees of phase sampling. 

o As part of our optimisation routines we recover a rotation pe¬ 
riod of 7.0 days. This is consistent with the measured chro¬ 
mospheric and coronal activity levels of the star (log(R)^^) = 
-4.38 + 0.05, this paper; log(Lx) = 29.02 ± 0.06, Kashyap 
et al. 2008). 

o The magnetic field reconstructions for HD 1237 are domi¬ 
nated by a band of strong uni-directional azimuthal field at 
high latitudes, accompanied by a complex multi-polar radial 
field distribution. The largest magnetic regions show field 
strengths of ~ 90 and 50 G for the azimuthal and the radial 
field components, respectively. 

o We note that the field topology recovered for HD 1237 is 
fully commensurate with studies of ^ Boo A and GJ 182, 
the two other stars that are closest to it in the mass, period 
diagram (Donati & Landstreet 2009). Larger sample sizes 
are needed to confirm these trends. If confirmed, it will be 
possible to predict the global magnetic field topologies and 
therefore the extended environments of planet-hosting stars 
at various stages of their evolution based on these fundamen¬ 
tal parameters. 


o We will address the influence of the different magnetic scales 
on the coronal structure and wind properties of HD 1237 
in a follow-up paper, using the ZDI maps presented here as 
boundary conditions (see Cohen et al. 2010, 2011a). Possi¬ 
ble star-planet interactions occurring in the system via tran¬ 
sient (e.g. CME events, Cohen et al. 201 lb) or quiescent phe¬ 
nomena (e.g planetary radio emission, Stevens 2005; Vidotto 
et al. 2012) can be also considered in future work. 

Acknowledgements. Based on observations made with ESO Telescopes at the La 
Silla Paranal Observatory under programme ID 089.D-0138. 
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Appendix A: Second Epoch Dataset (2012 Dec) 

We consider only the Milne-Eddington line profile for the 2012 
Dec dataset. Figure A.l contains the recovered ZDI maps for 
this epoch with the corresponding synthetic Stokes V profiles. 
Despite having a lower phase coverage in this case, we were 
able to recover robust magnetic field maps, fitting the spectro- 
polarimetric data up to an optimal reduced = 0.6 (see Sect. 
6.2). This small value of reduced results as a consequence of 
the fewer constraints available for this dataset. 



Fig. A.l. Results of the ZDI analysis for the second epoch observations 
of HD 1237 using the Milne-Eddington line profile. See caption of Fig. 
7 for more details. In this case the maps fit the spectro-polarimetric data 
to an optimal reduced = 0.6. In this case, ® = 0.0 is assigned to the 
observations acquired at HID 2456265.0. 


The field distribution clearly resembles the one obtained for the 
July dataset, with a large contribution from the azimuthal and 
radial components to the total field. The main large-scale mag¬ 
netic features are preserved between both observed epochs. This 
is consistent with the behaviour shown by the activity indica¬ 
tors (Sect. 4.2, Fig. 3) and the longitudinal magnetic field (Sect. 
5.1, Fig. 4) in the entire dataset. Some of the smaller-scale struc¬ 
ture is not recovered and the ring of azimuthal field is not clear 
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in this case. These changes are expected due to the number of 
observations and phase coverage in this epoch (e.g. Donati & 
Brown 1997). The slight shift in longitude is due to the initial 
phase selection in this case, where O = 0.0 is assigned to the 
observations acquired at HJD 2456265.0 (Table 2). 
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